Zero temperature metal-insulator transition in the 
infinite-dimensional Hubbard model 
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The Mott-Hubbard metal-insulator transition [hJJ3] is 
one of the most fascinating phenomena of strongly cor- 
related electron systems. This transition from a para- 
magnetic metal to a paramagnetic insulator is found in 
various transition metal oxides, such as V2O3 doped with 
Cr |3|]. The mechanism driving the Mott-Hubbard tran- 
sition is believed to be the local Coulomb repulsion U 
between electrons on a same lattice site, although the 
details of the transition should also be influenced by lat- 
tice degrees of freedom. Therefore, the simplest model to 
investigate the correlation driven metal-insulator transi- 
tion is the Hubbard model 10 pfl 
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The zero temperature transition from a paramagnetic metal to a paramagnetic insulator is investi- 
gated in the Dynamical Mean Field Theory for the Hubbard model. The self-energy of the effective 
impurity Anderson model (on which the Hubbard model is mapped) is calculated using Wilson's 
Numerical Renormalization Group method. Results for quasiparticle weight, spectral function and 
self-energy are discussed for Bethe and hypercubic lattice. In both cases, the metal-insulator tran- 
sition is found to occur via the vanishing of a quasiparticle resonance which appears to be isolated 
from the Hubbard bands. 
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various authors P,pR| p"2[. The criticism is partly based 
on the fact that the IPT is essentially second order per- 
turbation theory in U (although iterated due to the self- 
consistency appearing in the DMFT) whereas the metal- 
insulator transition happens at [/-values of the order of 
the bandwidth. 

Non-perturbative methods are clearly needed to clar- 
ify the situation. At finite temperatures, the Quantum 
Monte Carlo method (QMC) should give reliable results 
and recent QMC calculations by J. Schlipf et al. |T|] 
gave no indications for a first order transition at finite T. 
The experimentally found first order transition in certain 
transition metal oxides can therefore not be due to elec- 
tronic correlations as modeled in (1) and lattice degrees 
of freedom will certainly play a role at the transition. 

At zero temperature, the isolation of the quasiparti- 
cle peak and the appearance of a 'preformed gap' has 
been shown by S. Kehrein pT| to be in contradiction to 
a skeleton diagram expansion. Also, no 'preformed gap' 
has been seen in calculations based on the Random Dis- 
persion Approximation (RDA) J12| where the opening of 
the gap and the vanishing of the quasiparticle peak was 
found to happen at the same critical U. The results for 
the gap and the quasiparticle weight are obtained in the 
RDA from finite size scaling of Exact Diagonalization re- 
sults of clusters with up to 14 sites. 

Both QMC and RDA are non-perturbative approaches 
which can, in principle, be applied for arbitrary low tem- 
peratures (the resolution of the low-frequency behavior 
in QMC and RDA is limited by the number of time 
slices or number of sites, respectively). The only non- 
perturbative approach which is presently able to cover 
the very low temperature regime directly in the thermo- 
dynamic limit is the Numerical Renormalization Group 
Method (NRG). This method has been introduced by 
Wilson for the Kondo problem H| and has been ap- 
plied by Krishna-murthy et al. to the impurity Anderson 
model |15j . It has been later shown that the NRG allows 
for a very accurate calculation of dynamical properties 
of various impurity models J16|,pT[ . This is important be- 
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where c ia . [c% a ) denote creation (annihilation) operators 
for a fermion on site i, t is the hopping matrix element 
and the sum X)<i 7 > i s restricted to nearest neighbors. 
Despite its simple structure, the solution of this model 
turns out to be an extremely difficult many-body prob- 
lem. The situation is particularly complicated near the 
metal-insulator transition where U and the bandwidth 
are roughly of the same order and perturbative schemes 
(in U or t) are not applicable. 

With the recent development of the Dynamical Mean 
Field Theory (DMFT) ||-[| a very detailed analysis of 
the phase diagram of the infinite dimensional Hubbard 
model became possible. The Rerative Perturbation The- 
ory (IPT) results of || gave a first order metal-insulator 
transition at finite temperatures. The transition is occur- 
ing within a coexistence region of metallic and insulating 
solutions extending from T = up to T* sw 0.02IF (W: 
bandwidth). On approaching the metal-insulator tran- 
sition from the metallic side (i.e. on increasing U), the 
authors of H found a quasiparticle peak with vanish- 
ing spectral weight which becomes isolated from upper 
and lower Hubbard bands. A consequence of this result 
is that the opening of the gap and the vanishing of the 
quasiparticle peak do not happen at the same critical U. 
The possibility of this scenario has been questioned by 



cause in the DMFT the self-energy (or equivalently the 
single-particle spectral function) of an effective impurity 
model has to be calculated in the full frequency regime 
(for first applications to the Hubbard model see ]l8| , ^9| ). 
One therefore expects that the NRG gives equally ac- 
curate results for the effective impurity Anderson model 
appearing in the DMFT. However, due to the lack of ex- 
act results for, e.g., the metal-insulator transition in the 
Hubbard model, this cannot be proven so far. 

Here we concentrate on the Mott-Hubbard metal- 
insulator transition at zero temperature and half-filling. 
At T = 0, this transition is usually hidden by the tendency 
of the model to form an antiferromagnetic groundstate 
(as long as no frustration by, e.g., longer range hopping 
is included) . The results are discussed for both the Bethe 
lattice with infinite coordination number and the infinite 
dimensional hypercubic lattice. The hopping matrix ele- 
ment in the hamiltonian ([j]) is scaled as t — t*/y/Z with 
Z the number of nearest neighbors. In the following, we 
set t* = 1 as the unit for the energy scale. The resulting 
free densities of states for Bethe and hypercubic lattice 
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The effective bandwidth W = AJ j dep{e)e 2 is W — A for 

both p B and p^ c (the factor 4 is chosen so that the W 
corresponds to the actual bandwidth of the semi-elliptic 
density of states p&). 

Fig. [l]a shows the ^/-dependence of the quasiparticle 
weight 



Z = 



1 



1 _ 9SeS(t£) i : 

1 a. . \lj— 



(4) 



du 



for both lattices. Despite the different lattice structure, 
the critical value of U is approximately the same for both 
Bethe and hypercubic lattice, C/ Cj b ~ 5.88 = 1.47W and 
£4,hc ~ 5.80 = 1.45W. The different behavior of the 
Z(U)-curves for small values of U can be understood from 
second order perturbation theory which gives Z(U) = 
1.0 - 0.082 U 2 + 0(U 4 ) for the Bethe lattice jl| and 
Z(U) = 1.0 - 0.12 U 2 + 0{U 4 ) for the hypercubic lattice. 

Fig. |l|b shows the NRG-result for the Bethe lattice to- 
gether with results from calculations using the Random 
Dispersion Approximation (RDA) @ and the IPT ||. 
The NRG and RDA results agree very well up to £/«2.5 
but the NRG gives a long tail in the Z(/7)-curve end- 
ing at a critical U c which is considerably larger than the 
^c,rda ~ W. In my view, the difference for U > 2.5 
may be a consequence of the small system sizes presently 
taken into account in the RDA. 

The Z(U)-cuive from the IPT starts to deviate from 
the NRG result already for very small values of U. It 
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FIG. 1. U dependence of the quasiparticle weight Z; a: 
comparison of the NRG results for Bethe and hypercubic lat- 
tice and b: comparison of the NRG results with results from 
the RDA and the IPT (for the Bethe lattice). Also shown is 
the result for the critical U from the PSCM. 



has been found earlier that the critical U is overesti- 
mated by the IPT (C/ cJPT « 1.65W @) and that the 
Projective Self-consistent Method (PSCM) [g,|o| gives a 
lower C/ Cj ps ~ 1.46W which is in remarkable agreement 
with the NRG-result. The Z(C/)-curve from the QMC for 
small finite temperatures and a Bethe lattice density of 
states (not shown here) agrees well with the NRG result 
for T = up to U w 4.5 @. The critical values ob- 
tained from the QMC (e.g. U C:QM c = (1.26±0.01)W for 
T=l/30) are smaller than those from the NRG and the 
U Ct QMc(T)-curve shows a negative slope, a consequence 
of the higher spin-entropy of the insulating phase. 

The spectral functions A(uj) for Bethe and hypercubic 
lattice are compared in Fig. || for U — 0.8U C , U — 0.99U C 
and U—l.lUc (for details of the numerical calculations, 
see [|l9|). Although the semi-elliptic density of states pe 
is confined to the interval [—2,2], whereas the gaussian 
density of states phc has no cutoff, the structures ap- 
pearing in the spectral functions are very similar. In the 
metallic phase (for large enough values of U) the spec- 
tral function shows the typical three-peak structure with 
upper and lower Hubbard bands centered at ±U/2 and 
a quasiparticle peak at the Fermi level. For C/ = 0.99C/ C , 
the quasiparticle peak in both Bethe and hypercubic lat- 
tice seems to be isolated (within the numerical accuracy) 
from the upper and lower Hubbard bands, similar to what 
has been observed in the IPT calculations for the Bethe 
lattice ||. Consequently, the gap appears to open dis- 
continuously at the critical U. Note that, due to the 
broadening of the spectra JhJ , an accurate resolution of 
the high energy features, e.g., the band edges of the Hub- 
bard bands, is not possible. 
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FIG. 2. Spectral functions for Bethe and hypercubic lat- 
tice for various values of U. In both cases, a narrow quasipar- 
ticle peak develops at the Fermi level which vanishes at the 
critical U c - 



The term 'preformed gap' is frequently used to de- 
scribe the behavior seen in Fig. [2] although it is not clear 
whether a 'preformed gap' only means a strong suppres- 
sion of spectral weight between the Hubbard bands (as 
seen in both NRG and IPT) or an exact vanishing of 
the spectral function in a finite interval. Using the latter 
definition, it is not possible to decide within a numerical 
approach (like NRG, QMC or IPT) whether the system 
shows a 'preformed gap' or not. 

The three-peak structure in the spectral function and 
the isolation of the quasiparticle peak near the transition 
have important consequences for the behavior of the self- 
energy. Fig. shows the imaginary part of the self-energy 
for the same parameters as in Fig. 0. In the insulating 
regime (U = 1.1U C ) the self-energy has a pole at zero 
frequency £(z) = j + £ rcm (z) (z — uj + i0+, E rem de- 
notes the remaining part of the self-energy). There are 
several possibilities how the 1/z-term in £(z) develops 
when the transition is approached from the metallic side. 
One possibility would be that a is zero at the transition 
(U = U c ) and then increases continuously with increas- 
ing U. This imposes some constraints on the form of the 
spectral function because the weight a of the pole in the 
insulating regime is given by 
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If one assumes a powerlaw A(uS) oc \uj\ r for small uj, the 
exponent has to satisfy r < 1 for a to vanish at the 
transition. 

The other possibility is what is seen in the NRG-results 
for the Bethe and hypercubic lattice in Fig. g. The 1/z- 
term emerges from a two-peak structure in the imaginary 
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FIG. 3. Imaginary part of the self-energy for Bethe and 
hypercubic lattice for the same parameters as in Fig. ti. 



part of the self-energy. The weight of the peaks is roughly 
independent of U, while the position and the width van- 
ish with U — > U c (the position is proportional to \fZ s ). 
At the transition, the two peaks collapse and give rise 
to a single pole with weight a (these features have al- 
ready been discussed in |ll| , pi[ ; note that the two-peak 
structure in ImE does not imply the existence of poles 
in £(z)). This behavior is common to both lattice types 
studied here and it is only the [/-dependence of the width 
x, that differs between Bethe and hypercubic lattice. 

Note that the vanishing of the quasiparticle peak in the 
standard single impurity Anderson model (which occurs 
for U/ttA — > 00 [g2[) is also associated with the collapse 
of a two-peak structure in the self-energy. This is ob- 
served both in the wide-band and the narrow-band limit 
(see [E3[ for a discussion of the latter case) . 

The two-peak structure is related to the typical three- 
peak structure (quasiparticle peak plus upper and lower 
Hubbard bands) in the spectral functions for both single 
impurity Anderson model and the infinite dimensional 
Hubbard model. In both models, one has the relation 
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with the hybridization function A(z). The self-energy 
develops peaks at the frequencies where real- and imagi- 
nary part of G(z) are small, which is the region between 
the quasiparticle peak and the Hubbard bands. There- 
fore, all calculations for the infinite dimensional Hubbard 
model which give a well-pronounced three-peak structure 
in A(oj) necessarily produce the two-peak structure in the 
self-energy (examples are calculations from the QMC B 
and the Non Crossing Approximation |24[ ] ; at finite tem- 
peratures, the two-peak structure is broadened). 

We now turn to an additional feature seen in both IPT 



and NRG calculations: the coexistence of metallic and in- 
sulating solutions in an interval U Ct i < U < U c .2- Start- 
ing from {7 = 0, the metal to insulator transition occurs 
at the critical U c ,2 with the vanishing of the quasiparti- 
cle peak. Starting from the insulating side, the insulator 
to metal transition happens at U Ci i < U c p (the NRG 
and IPT give U c ,i ~ 1.25W for the Bethe lattice and the 
NRG gives U c ,i sw 1.15W for the hypercubic lattice). 

The coexistence of metallic and insulating solutions is 
probably connected to the structure of the self-energy at 
T = 0. When U is reduced below U c< 2, the <5-function peak 
in the imaginary part of the self-energy does not split into 
the two-peak structure which is found for the metallic 
solution. The <5-function peak only vanishes when its 
weight 7ra vanishes which happens at a lower value of U 
(the U c ,i). The [/-dependence of a near U c ,i is difficult to 
determine and it is presently not clear whether a vanishes 
continuously or not. 

The physical solution of the DMFT equations in the 
coexistence region is the one with the lower energy which 
turns out to be the metallic one, in agreement with |9|] 
(near the transition, the energy difference becomes too 
small to decide which solution has the lower energy, so 
that a small uncertainty remains near U c ^)- The coexis- 
tence of solutions therefore does not play a role at T = 
and can be neglected. 

In conclusion, we have investigated the zero- 
temperature metal-insulator transition in the Hubbard 
model for both Bethe and hypercubic lattice using a non- 
perturbative approach, the Numerical Renormalization 
Group method. The NRG calculations show that the 
details of the transition are very similar in both cases de- 
spite the different lattice structure. In the Bethe lattice 
case, the result for the critical U is in remarkable agree- 
ment with the result from the Projective Self-consistent 
Method. The NRG-results (in particular, the two-peak 
structure in the imaginary part of the self-energy near 
the critical U) cannot be explained within a skeleton di- 
agram expansion as shown in flllfl . Potential problems 
due to the fact that the derivation of the DMFT is based 
on such a skeleton diagram expansion (23| have still to 
be clarified. 

In order to bridge the gap between T = and the 
lowest temperatures accessible to the QMC-method, the 
NRG has to be extended to finite temperatures (work on 
this is in progress). This will allow to study the metal- 
insulator transition in the whole temperature range using 
non-perturbative methods. 
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